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We consider the directed network (DN) of edge states on 
the surface of a cylinder of length L and circumference C. By 
mapping it to a ferromagnetic superspin chain, and using a 
scaling analysis, we show its equivalence to a one-dimensional 
supersymmetric nonlinear sigma model in the scaling limit, 
for any value of the ratio L/C, except for short systems 
where L is less than of order C^^'^ . For the sigma model, 
the universal crossover functions for the conductance and its 
variance have been determined previously. We also show 
that the DN model can be mapped directly onto the ran- 
dom matrix (Fokker-Planck) approach to disordered quasi- 
one- dimensional wires, which implies that the entire distribu- 
tion of the conductance is the same as in the latter system, 
for any value of L/C in the same scaling limit. The results of 
Chalker and Dohmen are explained quantitatively. 
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I. INTRODUCTION 

Disordered conductors have been at the focus of ex- 
perimental and theoretical research for quite some time. 
Even properties of a single electron in random potential 
are quite nontrivial. One of the challenging problems 
still open in this area is the description of the transition 
between the plateaus in the integer quantum Hall (QH) 
effect. Chalker and Coddington introduced a network 
model to deal with this problem, and studied it numer- 
ically. Later several authors mapped this model to an 
antiferromagnetic spin chain, using replicas or supersym- 
metry to average over the disorder . The analysis of 
this spin chain is still complicated, and this model is not 
solved at present. 

Recently Chalker and Dohmen |^ introduced a three- 
dimensional (3D) version of the network model to study 
the possibility of QH effect in 3D conductors. This 3D 
network models a conductor composed of stacked coupled 
layers perpendicular to a strong magnetic field. Each 
layer may separately exhibit the QH transition and is rep- 
resented by the usual 2D network. In this system there 
exists a phase in which each individual layer is in the 
middle of a QH plateau, and all the electronic states in 
the bulk are localized. The only current-carrying states 
are the edge states on the surface of the conductor. These 



chiral surface states were studied in |p|-[10| . In the frame- 
work of the network model, these surface states form a 
2D directed network (DN), equivalent to one studied be- 
fore ||ll| . A distinguishing feature of the model studied 
by Chalker and Dohmen is the periodic boundary con- 
dition on the edge of each layer, making the system the 
surface of a cylinder. 

In this paper, we study (following Chalker and Doh- 
men) the conductance properties along the axis of the 
cylinder in the DN model. The paper is organized as fol- 
lows. We set up a general formalism for the description of 
the DN model using supersymmetry and the approach of 
Ref. II in section ||. Then we go to a continuum limit. 
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and map the DN model to a spin chain in section 
(this was also done in Refs. |,0,|l). Unlike the case of 
the QH transition, this spin chain is ferromagnetic. This 
allows us to analyze its properties using ideas of a recent 
paper on continuum quantum fcrromagnets jl^ . Scaling 
arguments, similar to those of Ref. |l^, show that all 
the moments of the conductance and other observables 
behave universally in the scaling limit (i.e. are given by 
universal scaling functions of dimcnsionlcss combinations 
of the hare couplings of the continuum model). 

In particular, we show that when the ferromagnetic 
chain under consideration is in the classical regime (called 
"renormalized" classical in Ref. |l^], but we drop the 
modifier because, due to the no-scale-factor universal- 
ity p^ , no coupling constants are renormalized), it can 
be further reduced to a ID classical non-linear sigma 
model, studied before by Mirlin et al. in the context of lo- 
calization in quasi-lD wires |ll3|] . Based on this reduction 
we show that in this classical scaling regime all the mo- 
ments of the conductance g are the same for both models. 
Thus, borrowing results of Ref. we can fully deter- 
mine the crossover functions for the mean and variance 
of the conductance of the DN model for any ratio of the 
length and circumference of the cylinder in the scaling 
regime. This classical scaling regime does not include 
the regime of very small L < 0{C^^^) that was termed 
OD in Ref. g. 

A recent paper [ p^ shows how one can in principle ob- 
tain the full probability distribution of the transmission 
eigenvalues from the non-linear sigma model, and proves 
that this distribution is identical to the one obtained from 
the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation 
| p5| , p^ of the Fokker-Planck (FP) approach which de- 
scribes the universal behavior of localization in quasi- 
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ID wires. In view of this equivalence, the probabiHty 
distribution of the conductance of the DN model in the 
classical scaling regime is the same as that of the quasi- 
ID model. For the quasi-lD model, a lot is known ex- 
actly about this distribution , and thus a nearly com- 
plete description of the conductance properties of the DN 
model is available. In section IV we give a direct argu- 



ment that shows how the DN model is related to a quasi- 
lD model, and thus to the DMPK equation. We compare 
our results with those of Chalker and Dohmen |^. We 
conclude in section 



II. THE DIRECTED NETWORK MODEL: 
GENERAL SETUP AND THE SYMMETRY 
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The DN is shown in Fig. |i] and consists of links and 
nodes. The links carry complex fluxes, and the nodes 
represent (unitary) scattering matrices S connecting in- 
coming (i,?') and outgoing (o, o') fluxes: 



= S 
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The scattering amplitudes a, . . . ,S correspond to elemen- 
tary scattering events shown on the right in Fig. For 
the time being they are assumed to be arbitrary complex 
numbers different for different nodes, which allows us to 
formulate our model for disordered samples with any re- 
alization of disorder. A particular distribution for the 
scattering amplitudes will be specified later. The verti- 
cal direction in Fig. is along the circumference of the 
cylinder. Later it will play the role of imaginary time for 
the spin chain, so we call the vertical coordinate t. In 
the T-direction the network has the size C = NtCLt where 
Nt is the number of "channels" through the system and 
ttr is a microscopic scale of the order of the mean free 
path of electrons. We impose periodic boundary condi- 
tions in this direction. In the x-dircction the network 
has finite length L = Na^, where N is the number of 
layers (or "sites" ) , and ax is the distance between them. 
In this direction the system is connected to ideal leads. 
In Fig. the edge states are numbered from 1 to iV = 5, 
and Nr = 3. 

Different correlation functions may be defined for this 
model and each of them may be represented in the first 
or second quantized way. As an illustrative and impor- 
tant example, we derive expressions for the conductance. 
The dimensionless conductance is given by the Landauer 
formula g — ix t^t = J2i f l^i/l^j where t is the total 
transmission matrix (with matrix elements tif) between 
left and right boundaries of the system. In the first quan- 
tized language, tif is given by the sum over "retarded" 
paths connecting an incoming link at on the left bound- 
ary and an outgoing link at Tf on the right boundary. 
Each path follows links only in the direction of the ar- 
rows and its contribution is the product of the scattering 



FIG. 1. The directed network (DN) model. The particles 
move on the links in the directions shown by the arrows. The 
bold line represents a typical path. The scattering amplitudes 
on such a path are as shown on the right. 

amplitudes along the path. One such path is shown on 
Fig. with bold lines. Similarly, t*j: is given by the sum 
over "advanced" paths where each node contributes the 
complex conjugate scattering amplitude. 

In the second quantized language, the sum over paths 
may be written as a trace of an "evolution" operator in 
a Fock space of bosons and fermions. To represent Uf 
{t*f) we introduce a retarded (advanced) boson (hi) 
and fermion fi [gi) on each link i. The numbers of bits 
of paths on each link play the role of occupation numbers 
of these bosons and fermions, and the collection of these 
numbers on horizontal section at a given r-coordinate 
specifies a state in the Fock space. Then the evolution 
operator V12 for a single node between the sites 1 and 2, 
which evolves quantum states on these sites by one step in 
the T-direction, is given by the sum of the contributions 
of all possible scattering events, described by iS. In a 
typical event involving only retarded bosons, k bosons are 
transferred from site 1 to site 2 with amplitude 7, 1 bosons 
are transferred from site 2 to site 1 with amplitude /3, and 
the remaining bosons stay on their respective sites, each 
contributing factors a for site 1 and 5 for site 2. This 
event gives the term {-f^l3^/k\l\){alf{a\Ya'^-'^5'^'^-'a\a\ 
(where Uai = a\ai etc.) in the evolution operator Vi2. 
After some rearrangement, the sum of all such terms for 
all the bosons and fermions may be written as 

^12 = : expf ^(4ai + flh) + ^{a\a2 + flh) 
\a 



+ J:(&2^i + 32.9i) + ^(^1^2 + 5132) j : 

where colons stand for normal ordering. 

The contribution of the boundary nodes connected to 
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the leads is different. When representing every 
leftmost node at r 7^ has only paths reflected off the 
left boundary (with the corresponding amplitude 6 or 5*), 
because only such paths contribute to Then the 

only possible event at such a node is that all the particles 
stay on the site 1, each contributing factors S for retarded 
and S* for advanced ones. The corresponding evolution 
operator is simply 

For the leftmost node at we also need to inject one 
retarded and one advanced path into the system. This 
is represented by the event where we create additional 
retarded and advanced particles on the site 1. For defi- 
niteness we choose them to be bosons. The correspond- 
ing evolution operator is |7pa|6|(S""i+"^'i = 
|7paj6|Vbi. Similarly, for the rightmost nodes at r ^ r/ 
we have 

Vn,n+i = 0"""+"^" , (4) 

and for the boundary node at Tf the evolution operator 
is \j\'^VN,N+iaNbN- 

The total evolution operator U is composed of all the 
Vi^i+i's in the following manner. Assume for definiteness 
that the number of sites N is odd. Then in our system we 
have even rows of links at integer times 0,ar,2ar, etc., 
where fluxes enter the system from the left and exit it 
to the right, and odd rows of links at half- integer times 
aT-/2, 3aT-/2, etc., where fluxes enter from the right and 
exit to the left. For even rows we form the product 

Veir) = VoiiT)V23iT) . . . VW-i,w(r), (5) 
and for odd ones 

Vo{t) - Vu{T)V3i{T) . . . VN.N+iir). (6) 
The operator U is then given by the product 

Nr-l 

U= Yl Vo[{n+l/2)ar)Ve(nar), (7) 

which is ordered with the earliest times at the right. Note 
that the only r-dependence in the operators Ve, Vo is 
through the iS-matrix elements, which implicitly depend 
on i and r. 

With the help of the operator J7, the conductance is 
given by 

g = STr^r^ \lr?ci\{n)h\{n)\lf?ciM{Tf)hN{Tf) u\ 

Ti,Tf 

(8) 

Here r.; and Tf label the times at which the creation 
and annihilation operators act on the states, 7^ and 7/ 



are the S matrix elements at the corresponding nodes, 
and T orders times, placing the earliest at the right. 
"STr" stands for the supertrace in the Fock space, which 
weights all the states with the factor of (—1)^^, where 
Np = X]i3=i("'/» ~^^gi) is the total number of fermions in 
a state, STr(. . .) = Tr((-1)^^ . . .). That is, states with 
an odd number of fermions contribute to the sum with 
a negative sign. This, together with the periodicity in 
the r-direction, ensures the cancellation of bosonic and 
fermionic contributions from closed paths not connected 
to the leads. These closed paths were, as usual, absent 
from the original first quantized formulation of the prob- 
lem. 

In Eq. (^), Ti and r/ take only half-integer values ar/2, 
Sar/^, so on, because we create and destroy bosons only 
on half of the links belonging to the sites 1 and N. Using 
the commutation relations between the bosonic operators 
and 1^'s, we can rewrite the expression for g in the form 

g = STr(TY,era\{r)b\{r)Y,er'aN{r')bN{T')u\ 

(9) 

where now t and t' take all possible values, but = 1 
(— 1) for integer (half-integer) values of r/ar- 

This is a good place in which to discuss in more detail 
the choice of boundary conditions for our model. In the 
x-direction the presence of the ideal absorbing leads at 
the boundaries of the system means in the first quantized 
language that we do not include contributions from the 
paths leaving or entering the system unless we specifi- 
cally calculate some correlators between the boundaries 
(like conductance g). In the second quantized language 
this translates to the following. We can imagine having 
two additional vertical sets of links in the leads (0th and 
N + 1st sites) on which we have no bosons or fermions, 
so these sites always carry the vacuum state |0). With 
this constraint the boundary operators Vbi and Vn.n+i 
are seen to be special forms of general Vi^i+i^s, Eq. (||), 
acting on the vacuum at the left or right. Later this con- 
straint on the states at the boundaries will give, in the 
T-continuum formulation, the symmetry-breaking term 
in the Hamiltonian and in the action, and will fix the 
boundary conditions in the continuum field theory, see 
Eq. d^) below. 

We now discuss the supersymmetry properties of our 
formulation. For each site i we can form 16 bilinears 
in our bosonic and fermionic operators. They represent 
the 16 generators of the Lie superalgebra u{l, 1 \ 2). We 
arrange these generators in a 4 x 4 matrix, or "super- 
spin" Ji'. 
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J,. = 



/ t 1 /r. 

'aia] - 1/2 




aibi 




M 


hfl - 1/2 


fibi 








-b\br - 1/2 


-by 


\ 9}al 






ajg^ - 1 



1/2/ 
(10) 



We adopted "advanced-retarded" arrangement of the 
generators rather than the "boson-fermion" one (see ^] 
for details) . In this scheme the diagonal 2x2 blocks of Ji 
contain the generators of the subalgebra u{l \ 1)®m(1 | 1). 
Then one can show that the operator V12 commutes with 
all 16 components of Ji + J2 {and thus has w(l, 1 1 2) as 
the symmetry algebra), 



iJl+J2)Vi2^Vi2{Jl+J2), 



(11) 



if and only if the scattering matrix S is unitary: S^^ = 
5^ (the analogous result for the Chalker-Coddington 
model was found in Ref. [^). Thus, the products of 
the l^i^i+i's in Ve and Vo commute with the total su- 
perspin X^iLi Ji for any realization of disorder, provided 
that all the scattering matrices are unitary, except for 
the boundary operators Vbi and Vn,n+i- These bound- 
ary operators commute only with generators from diag- 
onal blocks of Ji and Jat, correspondingly, and break 
the symmetry of the total evolution operator U down to 
u{l 1 1) © u{l 1 1). Therefore, for any realization of the 
disorder, unitarity of the scattering matrices ensures the 
global u(l, 1 1 2) symmetry of our problem, which is bro- 
ken only by the boundary constraint. The use of the su- 
pertrace, not the ordinary trace, is essential in maintain- 
ing the supersymmetry in the presence of the periodicity 
in the r-direction. 

In terms of the superspins, Eq. for the conductance 
may be rewritten as 
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STr( Tj2srJ!\r)J2eA-irJhHr')u). (12 



Here Jj^ = —b[a[,Jj^ = a^bN are particular compo- 
nents of superspins Ji and Jn at the boundaries of the 
system. The alternating sums of the type 



X. = (-l)'^£.jf(r) 



(13) 



represent the total current through the system in the x- 
direction. Using the supersymmetry of the operators V, 
Eq. (pr]), we can show that this current is conserved, i.e., 
correlators of X,;'s do not depend on the position labels i, 
except at the boundaries. We see then that Eq. ( p^ is 
the usual Kubo-type formula, relating the conductance 
to a current-current correlator, g = — STr(TIiX,[/), for 

Similarly we can write down expressions for the mo- 
ments of the conductance or other correlators. For mo- 
ments higher than the second, we need to introduce some 



additional structure. Namely, we introduce n replicas of 
our bosons and fermions and sum the bilinear operators 
in the exponentials in the expressions for Vi,i+i's over 
replica indices. This is necessary, because to represent 
the TO-th power of the conductance we need to create 
at the boundary m different retarded and m different 
advanced particles, and we need n > to/ 2. The corre- 
sponding model has global u{n,n \ 2n) symmetry broken 
down to u{n \ n) © u{n \ n) by the boundary constraint. 
We should point out that in this case replicas are not 
introduced to average over the disorder, and their num- 
ber n is not taken to zero in the end. In this model, 
results for the m-th moment are independent of n, pro- 
vided TO < 2n, because any "excess" replicas cancel by 
supersymmetry. For simplicity we continue to assume 
n = 1. 

So far we did not have to specify the nature of the 
disorder. Now we assume a particular distribution of the 
scattering matrices, used before in 0,^. Namely, we take 
every scattering matrix to consist of a product of two 
diagonal unitary matrices with a real orthogonal matrix, 
with entries 



/3 = 7 = t>0, 



5={l-t') 



2N1/2 



(14) 



in between them, the latter matrix being the same for 
all nodes. The phases from the diagonal unitary ma- 
trices are associated with links rather than with nodes 
and are assumed to be uniformly distributed between 
and 27r, independently for each link ||l|,^. Thus averaging 
over disorder corresponds to integration over all the link 
phases. 

It is easy to see that this averaging produces a local 
constraint on the states in the Fock space. Let each 
retarded bit on a given link at the i-th site at time t 
contribute the factor e*''', and each advanced one the 
factor e"*"^, to the amplitude associated with a given 
path. Then if we have to retarded and n advanced 
bits on the link, averaging over random phase Lp gives 
exp{i{ra — n)Lp)d(p — Sm,n, that is, we get zero unless 
the number of retarded and advanced bits is the same. 
Thus, we have a local constraint on the number of bits of 
paths on each link. In the second-quantized language this 
means that the averaging projects our evolution operator 
to the subspace specified by riai+Ufi = nu+ngi for each i 
and each time r. This subspace forms a highest weight 
irreducible representation of the algebra u(l, 1 1 2), with 
the vacuum |0) for the site i being the "highest weight" 
vector. This representation was first obtained for the 
Chalker-Coddington model [Q in [|| and was also dis- 
cussed in [Q. 
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III. MAPPING TO A SPIN CHAIN AND 
SCALING ANALYSIS 

So far in our derivation we retained discreteness of the 
DN. Unhke in the other pubHshed derivations (5|j^,P, we 
did not have to take the time-continuum Hmit from the 
start. However, we do so now in order to obtain a contin- 
uum field theory for our system and to study its universal 
scaling properties. In the r-continuum limit we assume 
i <C 1 and expand the Vi^i+i's, projected to the con- 
strained subspace, to second order in t. The result of the 
expansion is = 1 — t'^stiJiJi+i + C(f*), where Ji 

are superspins in the highest weight representation men- 
tioned above. Here "str" stands for the matrix supertrace 
in the space labeled by the upper indices of Jf*. In other 
words, for any 4x4 supermatrix A, str A = tr rjA where 
"tr" is the usual matrix trace and r/ is a diagonal matrix 
with entries (1, —1, 1, —1). Now we combine all the V's 
and reexponentiate. We also replace sums over r by inte- 
grals. The result is the evolution operator in imaginary 
time r of a ID quantum ferromagnetic spin chain with 
Hamiltonian 



e /^"^ 1 
H[J] = — 51 strJ, J,;+i + -strA(Ji + Jn) 



(15) 



Here A is the diagonal matrix diag(l, 1,-1,-1). The last 
term in H comes from the boundary nodes. It may be in- 
terpreted as saying that there are, at the boundaries, two 
additional spins Jq and Jn+i fixed to a particular "direc- 
tion" Jo — Jjv+i = A/2 by an infinitely strong magnetic 
field coupling via a Zeeman term in the superspin space. 
This is another manifestation of the boundary constraint 
mentioned above. 

Using the algebra of the generators in the same r- 
continuum limit, the expression (^3|) for the total current 
becomes 

^2 rC ^2 rC 

— / dr[Ji,Ji_i] = / dT[Ji,Ji+i]. (16) 

Jo o,T Jo 

For the particular components of the current at 
the boundaries these expressions reduce to Tf^ = 
(t2/2a^)/dT[Ji, A]3i = (t^/ar) JdrJf^ and, similarly, 
Ijf = ityar)JdTj],\ 

Next, we introduce supercoherent states and represent 
quantities of interest as path integrals over some super- 
manifold, see |^,|l^ for details. The resulting theory has 
the action 



S = Sf_ 



dTH[Q{T)], 



(17) 



where Sb is the Berry phase term (specified below), 
Qiir) is a supermatrix taking values in the coset space 
U{1, 1 1 2)/[/(l 1 1) X U{1 1 1) (where [/(1, 1 1 2) is the su- 
pergroup of which the Lie superalgebra is u(l, 1 1 2), etc.), 



and H[Q{t)] is obtained from H[J] by replacing every Ji 
with Qi{T)/2 (we make the same replacement in the ex- 
pressions for the total current). In the path integral, all 
the components of Q obey periodic boundary conditions. 
The difference, in the case of the fermionic components, 
from the usual path integrals for fermions, which obey 



antipcriodic boundary conditions |19|, is a direct conse- 
quence of the factor (— l)^'' in the definition of STr. 

The scaling properties of quantum ferromagnets were 
discussed in Following that paper, we take the spa- 
tial continuum limit of the action S and the current I^. 
The resulting continuum ferromagnet has the action 



L nC 

dx 
Jo 



du sir Q{u) 



dQ{u) dQju) 
du dr 



-str(VO)' 



and the current becomes 



I{x) = -2ps / drQVQ 



(18) 



(19) 



Here Mo = l/2ax and ps = t^ax/^ar are the magneti- 
zation per unit length, and the spin stiffness in the fer- 
romagnetic ground state respectively. Also, in the Berry 
phase term, the first term in ^cont, Q('u) = Q{x,t,u) 
is some smooth homotopy between Q{x, r, 0) = A and 
Q{x,T,l) — Q{x,t) (for details see |18 ). The circum- 
ference of the cylinder C plays the role of the inverse 
temperature. The last term in H[Q{t)], or the interpre- 
tation Jo ~ Jn+1 = A/2, forces the field Q{x, t) to take 
the boundary values 



g(0,r) ^Q(L,r) = A. 



(20) 



A boundary condition of this form to represent ideal ab- 
sorbing leads is usual in the non-linear sigma model for- 
mulation of the theory of localization (to which, we note, 
the model (|l8|) is noi equivalent). Given that the effective 
supersymmetric spin system is ferromagnetic, and that 
fiuctuations in the length of the spins can be neglected, 
the form of the action Scant is dictated by the supersym- 
metry, which holds universally for any short range form 
of disorder. 

The action ( p^ clearly shows an anisotropic scal- 
ing which reflects the difference between the x- and r- 
directions in the DN. In the original discrete version of 
the DN model, if we consider the local behavior without 
the periodic boundary condition, the mean conductivity 
behaves diffusively in the a:-direction, with r playing the 
role of time, with the diffusion constant 



D = 



(21) 



as shown in Refs. (However, even without the pe- 

riodic boundary condition, less trivial behavior is found 
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for higher moments (T^.) In the spin chain language this 
anisotropy is manifest in the usual quadratic dispersion 
relation of the spin waves. The action jl^), linearized 
near the ferromagnetic ground state, describes (on ana- 
lytic continuation to real time) spin waves with the dis- 
persion uj = 2psk'^ /Mq. The difference between the ratio 
2ps/Mo = t'^a^/ur and the diffusion constant D is due to 
the particular order of time- and space-continuum limits 
taken to obtain the action (|lj). The more general long 
wavelength behavior {kxCix ^ 1, kT-a-r <^ 1), consistent 
with Eq. (pl|), should be described by the same action 
Scant but with the more accurate expression for the spin 
stiffness, 

In general, the continuum description requires that C 3> 
Or, L ^ a^. 

Although S'cont is a highly non-trivial interacting quan- 
tum field theory, it is nevertheless possible to make some 
simple exact statements about it. First, its "ground 
state" (which dominates the functional integral in the 
limit of zero "temperature" , C = oo) is fully polar- 
ized and fluctuationless, so {Q)c'=oo — A. As a result, 
the parameter Mq, which measures the Berry phase due 
to adiabatic changes in the ground state polarization, 
cannot have any non-trivial zero temperature renormal- 
ization Q. Second, not only the ground state, but 
some low-lying excited states are also known exactly: 
the linearized single spin-wave states are in fact exact 
eigenstates of the full Hamiltonian, and their disper- 
sion u! = 2psk'^ /Mq has no corrections from the non- 
linearities. This implies that also has no non-trivial 
zero-temperature renormalization . We also mention, 
as an aside, that the mapping from the DN model to 
Scant can be generalized to higher spatial dimensions d. 
As a result, the critical dimension d — 2 (for the proper- 
ties at large C) pointed out in Ref. (l2j is related to the 
same fact obtained in Ref. |]lT| . 

The absence of any non-trivial renormalizations of 
Mq and ps, the absence of any additional relevant op- 
erators in 1 + 1 dimensions, and the related absence 
of ultraviolet divergences in physical quantities in the 
1 -|- 1-dimensional quantum field theory, are responsi- 
ble for the phenomenon of "no-scale-factor universal- 
ity" [ p^|j2^ . Among its implications is that scaling anal- 
ysis of ^cont reduces to the naive analysis of engineer- 
ing dimensions. The field Q{x,t), being subject to the 
constraint Q'^{x,t) = 1, is dimensionless. Then simple 
power counting gives the engineering dimensions of the 
couplings: dim7\fo = (Icngth)^^, dim ps = length/time. 
The no-scale-factor universality also implies that all the 
observables are functions of dimensionless combinations 
of the bare couplings A/q and pa and the large scales, 
which, for the conductance, are just L and C . For exam- 
ple, the mean conductance is given by a function 
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FIG. 2. Sketch of the crossovers as a function of 
the length L. The first crossover Q is at the scale 
L ^ Lo = {C pa/MoY^'^ where the low-lying energy level split- 
tings of the quantum ferromagnet ~ PsIMqI? are of order the 
"temperature" 1/C. The second crossover is near the local- 
ization length 5 = 8Cps. The conditions for the validity of the 
continuum theory require ^ 3> -£/o, and so the two crossover 
scales are well separated. The quantum ferromagnet behaves 
classically for all L ^ Lq (but with no restriction on L/Cps) 
and is then described by the one-dimensional supersymmet- 
ric non- linear sigma model Si (Eq. (^)), which is therefore 
valid across the crossover at L ~ Cpa. The Fokker-Planck ap- 
proach of Section ^ has the same regime of validity {L ^ Lq, 
L/Cps arbitrary). The 2D chiral metal regime was labelled 
OD in Ref. 1^, but we prefer the present terminology for rea- 
sons discussed below Eq. (^l[). 



(5) = '^9 



L 



L 



(Cp,/Mo)i/2'Cp. 



(23) 



Clearly, we could have chosen other combinations of the 
two arguments of the ^g, and the physics behind our 
particular choices will become clear as we proceed. The 
continuum theory Scant applies to any particular lattice 
quantum ferromagnet provided the lengths L, C are suf- 
ficiently large. Specifically jl2|, we require MqL ^ 1 and 
PsMqC » 1 (or L > a^; and C > a^(l - t'^)/t'^, as well 
as C ^ fli-). Comparing with the function <I>g(£i,£2), we 
see that the condition on C implies that £i ^ £2, and <I>g 
becomes universal in this regime, with no non-universal 
factors in either the scale of $g or its arguments. No- 
tice that the arguments ^1, £2 can take all positive value 
< ^1, ^2 < 00 while satisfying £1 ^ £2- We also empha- 
size that the condition for universality may be satisfied 
for any ratio L/C, i.e. for short as well as long cylinders. 

We first present our results for the case L ^ Lq = 
{Cps/Maf/^. The length Lq was identified in Ref. [p|, 
and we will discuss the regime L < Lq later (see Fig ^). 
For L ^ Lq, we have from (E3[) 



(24) 



where ^gi{£) = ^g{oo,£). Notice that the ratio L/C is 
still allowed to be arbitrary. The form (24) implies that 
any localization length ^ must obey 1^ oc Cps, which was 
the basic result of Chalker and Dohmen which we 
have now explained. One might have been tempted to 
use the incorrect naive argument, that distances in the 
x-direction (such as ^) should scale as the square root of 
distances in the r-direction (such as C). However, under 
such scaling, ps is not scale invariant, while Cps does 
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scale as a length. In the underlying continuum quantum 
ferromagnet, there are two dimensionful parameters ps 
and A/o, so this is not a scale-invariant system either, 
unlike, for example, the usual 2D non-linear sigma model 
at the QH critical point. 

It is possible to obtain scaling functions like <I>gi di- 
rectly from a certain field theory. They are current- 
current correlators in the path integral with the ID action 



^1 = / dxsiT{VQf, 

^ Jo 



(25) 



obtained from Scont by neglecting the r-dependence 
of Q |2^]. The contributions of all modes with a non- 
zero "frequency" along the r direction to the coefficient of 
(VQ)^ in Si are suppressed by powers of Lq/C and Lq/L. 
Eqn ( |25| ) is exactly the action of the ID non-linear sigma 
model studied by Mirlin et al. in Ref. ||l3|. In that pa- 
per, the authors used harmonic analysis on superspaces 
to diagonalize the transfer matrix of this ID model. As 
a result they were able to express the mean conductance 
{g) and its variance var t/ as rather simple integrals/sums 
over spectral parameters. For example, {g) is given by 



(.9) -2 



m>0 
odd 



d\ tanh ^ 

2 m/ 



■ exp 



(m2 + A2)L 



2^ 



where f is the localization length. 



(26) 



(27) 



(Note that we do not include any spin degeneracy, unlike 
Ref. [§).) 

In the ID metallic limit of short cylinders, Lq "C i <C 
Eq. ( p6| ) reduces to [|j 



c e 



X = 

2L L l-t^ ar' 



(28) 



while in the opposite ID localized limit of long cylinders, 
X » ^, one obtains 



(5)L»e - 2 



2L 



3/2 



exp 



2^ 



(29) 



An expression similar to Eq. (|2q ) can also be written for 
the variance varg, see Ref. [ |13[ . In the metallic regime 
Lq ^ L ^ S^, this yields the known result for universal 
conductance fluctuations in ID, varg = 1/15. 

We now turn to a discussion of the regime L < Lq 
(Fig H). Different behavior arises here because Q, in 
the ferromagnetic language, the "temperature" 1/C be- 
comes of order the low-lying level splittings, which are 
of order Ps/MqL^. In contrast, the classical regime we 
have discussed above is where the temperature is much 



greater than these quantum mechanical splittings. The 
crossover from the very short regime, L <^ Lq (termed 
OD in Ref Hi), to the ID metallic regime, Lo <C L < 
will be described by the L/Cps — > limit of scaling forms 
like (^^. The scaling forms can be expanded as a per- 
turbation series in increasing powers of L/Cps, times a 
universal function of L/Lq in each term. For (g), we can 
show that this takes the form 



(30) 



To leading order in L/Cps there is no dependence on 
L/La, and the result ( pSj) is valid in both the OD and 
the ID metallic regimes. The possible term of or- 
der {L/Cpsf' vanishes identically, consistent with the 
known result that the leading "weak localization" cor- 
rection to {g) vanishes in the quasi-lD metallic regime 
L/Lq — > cxd in the present (unitary) case. The next 
term in the expansion in L/Cps in ( pO| ) does have a 
non-trivial crossover at the scale Lq, described by the 
universal function T, which should approach the known 
quasi-lD metallic result —1/180 as L/Lq oo jisj. In 
contrast, for the variance of g we expect 



var g 



(31) 



where the universal crossover from the OD to the ID 
metallic regime is now evident in the leading term T, 
which must tend to 1/15 as L/Lq oo. A result of 
this type, for varg, has been claimed by Mathur Q and 
Yu who find var 3 oc Cps/MoL"^ in the OD hmit 

L <^ Lq. This result resembles that in isotropic systems, 
for example in 2D in the metallic regime (5) ^ 1, where 
var g oc W/L (where W is the width), as W/L 00 
(see the results in Ref. [^). Thus a better name for 
this regime in the DN model would be 2D chiral metal, 
and the crossover at Lq = {Cps/Mo)^/^ is from 2D to ID 
behaviour. 



IV. FOKKER-PLANCK APPROACH AND 
COMPARISON TO EARLIER WORK 

The ID non-linear sigma model considered in the pre- 
vious section is well suited for obtaining moments of the 
conductance, but not for its distribution. However, as 
we pointed out in the introduction and as we now dis- 
cuss, this ID model is completely equivalent to the FP 
approach to the conductance of quasi-lD wires in the 
limit of infinite number of channels ("thick wire limit"), 
see Ref. [Q. Then wc can use the results of this ap- 
proach to obtain further properties of the DN model. 
In the FP approach one concentrates on the eigenvalues 
T„ = (coshxn)"^ of the transmission matrix t^t. The 
probability distribution of the parameters Xn satisfies 
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the DMPK equation, which was solved approximately 
in the localized and metallic regimes (see, for example, 
Ref. [^) and exactly (for the unitary case) in Ref. [ p7| . 
We summarize some of the results of this solution. 

In the metallic limit, L <§; ^, the Xn have statistical 
fluctuations, but the mean density of Xfi S IS uniform, 
and the mean conductance is equal to this density. This 
may be rephrased by saying that the parameters are 
equally spaced in the average, (a;,i)L<5 — n{xi)L<^(^, and 
the mean conductance is equal to the inverse of the first 
parameter (g)L<j = (xi)£;^^ = 

In the opposite localized limit, L ^ S^, the parameters 
Xji are self-averaging (with normally distributed fluctua- 
tions), and their mean positions are again equally spaced, 
but they are offset from the origin by a half of the spacing 
between them: {xn)L^^ = {n— 1/2) (2L/^). In the local- 
ized limit, the conductance is dominated by the smallest 
parameter, g « 4exp(— 2xi), and therefore its logarithm 
is normally distributed with 

2L 

{\ng)L»i^-2{xi)L»i^--^. (32) 

The product of (g)L<^ and (In g)L>^ is a universal num- 
ber: 

(ln5)L»5 = -1, (33) 

which is characteristic of the universal crossover from 
metallic to localized behavior. If we now introduce 
rescaled parameters as in i^n — (C/L)x„, then 
('^i)l>C = C/S., which for a given value of C differs by a 
factor of 1/2 from (i^i)L<{ = 2C/^. 

The equivalence, mentioned above, of the ID sigma 
model and the FP approach motivated us to look for 
a direct mapping from the DN model to the FP equa- 
tion, and thus for a better understanding of why the DN 
model behaves as a quasi-lD conductor. For this pur- 
pose instead of considering the evolution operators Ve 
and Vo , which are the transfer matrices in the r-direction 
(or "row transfer matrices" ) in the second quantized lan- 
guage, we should concentrate on the transfer matrices 
in the x-direction, or "column transfer matrices" . For a 
single node such a transfer matrix M. connects fluxes on 
the left {i, o) and on the right («', o') of the node and is 
given by 

with a, . . . , (5 being the same as in the Eq. (|]). After 
going from one site to the next through one column of 
nodes the directions of the fluxes are reversed at each r- 
coordinate. Then the natural transfer matrices Mi-i^i^i 
(which can be simply multiplied) are composed of all the 
A^'s for two adjacent columns of the nodes to the left and 
right of the i-th edge state. The total transfer matrix for 



the DN with an odd number of edge states A'' is then 
given by = M^-i^n+i ■ ■ ■ M2,4Mo,2- This matrix 
connects A^^-dimensional vectors of fluxes on the right of 
the system (/', O') with the ones on the left (/, O): 

(?)-'~(^)^('-»:)(o)- <») 

Now we neglect all the link phases (we will reinstate 
them later) and consider the limiting case a = S = 
and 7 = — /3 = 1, where back scattering from right- to 
left-moving flux, or vice versa, is absent. In this limit, 
the column transfer matrix Mi_i^i+i becomes the "shift" 
matrix, which means that all the right- (left-) moving 
fluxes are transferred without any change along the a;- 
direction by 2ax, and along the r-direction by Ot (— Ot)- 
In other words, Mi-i^i+i cyclically shifts all the right- 
moving fluxes by ar and all the left-moving fluxes by 
—a-r- It is now easy to see that, due to the periodicity in 
the T-direction, if we multiply together Nj- such column 
transfer matrices, starting with Mq 2i we get the identity 
matrix: M2Nr-i — 1- This is exactly the situation shown 
on Fig. 1^, where we chose N — — 1 = 5. 

When we introduce a small amount of back scattering, 
taking a — b <^ \ and 7 = — /? = 1 -|- 0{p?)^ the total 
transfer matrix M2n^-\ is still close to the identity ma- 
trix. However, the back scattering will produce nonzero 
off-diagonal elements. The paths contributing to the el- 
ements of blocks 7712 and ma of the matrix M2n^-\ will 
have at least one back scattering event on them, and, 
therefore, all the matrix elements of TO2 and TO3 will be 
of the order 0{a). Similarly, the off-diagonal elements in 
TOi and 7774 will come from the paths having at least two 
back scattering events, and will be of the order 0{o?). 

When we reinstate the link phases, the matrix ele- 
ments of the matrix M2n^ - 1 will also acquire some phase 
factors. These phases are correlated, because the paths 
contributing to two different matrix elements may have 
links in common. However, the resulting transfer ma- 
trix must be pseudo-unitary due to current conservation 
and can be factorized into a product, consisting of a real 
matrix sandwiched between two block-diagonal unitary 
matrices The real factor in this decomposition is of the 
form described in the previous paragraph. Thus, the 
back scattering is of the same order between any right- 
and left-moving channel on a mesoscopic scale 2Nrax, 
which is much smaller than the localization length ^ of 
Eq. (p7|), because in the limit we are discussing now 
t^/(l — t^) — (/3/a)^ » 1. This equal mixing among 
the channels, happening in the metallic regime L <g; ^, 
well before localization sets in, is the essential property of 
quasi-lD systems. By contrast, in isotropic 2D systems, 
like the original Chalker-Coddington model [Q, the back 
scattering mixes the channels only locally in the trans- 
verse direction, and when the system is at the QH critical 
point, the value of the length L at which the channels are 
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all equally mixed is of the order of the width W. For such 
systems with L ^ W localization sets in with localization 
length ^ ~ W, so there is no metallic quasi-lD regime. 
Such a regime occurs when a separation of scales ^ ^ W 
takes place, but only if the dimensionless 2D conductiv- 
ity (Jxx is large, whereas at the QH critical point axx 
takes on a universal value of order 1. For the DN model 
the crossover to ID behaviour occurs at L ~ Lq ~ C^/^, 
which is much less than the localization length ^ ~ C, so 
a ID metallic regime occurs. 

In quasi-lD wires (see for a review of the scatter- 
ing approach to disordered conductors), one starts with 
the idea that a wire of length L can be obtained by com- 
bining many short (mesoscopic) building blocks of length 
SL. The transfer matrix Ml for the wire is a product of 
the transfer matrices Msl for the individual blocks. The 
parameters Xn, introduced above, are simply related to 
the eigenvalues of the transfer matrices. Upon multipli- 
cation of the transfer matrices of the building blocks, the 
parameters Xn perform a random walk, and their proba- 
bility distribution for the total transfer matrix is obtained 
by repeated convolution of the distribution of x„ for the 
building blocks. Upon taking a limit, in which Msl — * 1 
as (5L 0, the distribution of the Xn is described by 
a universal FP equation with L a continuous variable, 
this is the DMPK equation. In our system, the building 
block is the cylindrical DN of length 2Nrax- Therefore, 
the universal properties of the DN model in the quasi-lD 
scaling limit should coincide with those of the DMPK 
equation, at least when t^/(l — i^) S> 1. 

Finally, we want to compare our results with those of 
Chalker and Dohmen Q . They introduced an amplitude 
ratio Aqd = ^cd / C and found it was equal to — t^) 
in agreement with numerical simulations, while from our 
Eq. (HI) we find for this ratio A = 2t^ax/{l - i^)ar. 
First, we should point out that Chalker and Dohmen 
measured their L and C in units of ax and a^., respec- 
tively. This takes care of the factor axja-r. We attribute 
the remaining factor of 2 difference between A and ^cd 
to the different conventions used in the definitions of the 
localization length i^. It appears from the equations of 
Chalker and Dohmen that they defined their localization 
length ^cD through the decay of the typical conductance 
in the localized regime, as (Iny) = — L/^cd, while we 
used the more conventional definition, such that Eq. (^) 
holds. Then our ^ — 2^cd, which explains the factor of 2. 

There is one more discrepancy between our results and 
those of Chalker and Dohmen. From their expressions 
for the conductance one finds that the universal num- 
ber of Eq. ( ^ ) is —2. This may be explained as follows. 
Chalker and Dohmen used rescaled parameters , intro- 
duced above, and made a conjecture (seemingly based 
on the ID model) that the value of v\ is the same in 
the metallic and localized limits. For the ID model this 
conjecture does not hold, as we saw above. If instead of 
this conjecture Chalker and Dohmen had assumed that 



the t'n's behave in the same way as in the ID model, i.e. 
that (j^i)l<5 = 2(i'i)l>{, as predicted by our mapping, 
their expressions would completely agree with Eq. (|3^). 
We conclude that our results agree with the numerics of 
Chalker and Dohmen, and explain why their heuristic 
argument works (after factors of 2 are corrected). 

V. CONCLUSION 

In conclusion, we have considered the directed network 
(DN) of edge states on the surface of a cylinder. We 
claim that, in a certain scaling limit, the DN is equiv- 
alent to the ID supersymmetric non-linear sigma model 
and to the random matrix model used before to describe 
the transport properties of quasi- ID wires. Using the 
known results for this ID model, we obtain a description 
of the conductance properties of the DN model in this 
scaling limit. In particular, we give exact expressions for 
the mean conductance (17) and the correlation length ^ 
of the DN model in the scaling limit for any value of 
the ratio C jL of the circumference and the length of the 
cylinder, Eqs. (^6|, |2^), while expressions for the variance 
and the distribution of the conductance may be found 
in the literature on quasi-lD wires ]l3| , [r7| , p3t . These re- 
sults are universal; in particular, they should not depend 
on the precise distribution chosen for the disorder (pro- 
vided its correlations are short-ranged). They are valid 
except for short systems with L Lq = ip PsIMq)^!"^ or 
less, which are expected to behave as 2D chiral metals. 
In essence, the crossover from 2D to quasi-lD behaviour 
is governed by anisotropic scaling, and hence occurs at 
L ^ Lq ^ (71/2^ while localization sets in at larger scales, 
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